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We derived analytic formulae for the polarization single-spin asymmetries (SSA) in the Bethe- 
Heitler process e~ + p — > e~ + 7 + p. The asymmetries arise due to one-loop QED radiative 
corrections to the leptonic part of the interaction and present a systematic correction for the studies 
of virtual Compton Scattering on a proton through interference with the Bethe-Heitler amplitude. 
Considered are SSA with either longitudinally polarized electron beam or a polarized proton target. 
The computed effect appears to be small, not exceeding 0.1 per cent for kinematics of current virtual 
Compton scattering experiments. 

I. INTRODUCTION 

Experiments on virtual Compton scattering (VCS) and deeply- virtual Compton scattering (DVCS) are an important 
part of nucleon structure studies at major electron-scattering laboratories, and at Jefferson Lab in particular [l|. The 
reason for such a high interest is that VCS allows to access 3-dimensional parton distributions of a nucleon and has 
a potential of resolving the role of parton orbital angular momentum in the nucleon spin problem, see the original 
papers and recent reviews 0, Q . 

In experimental observations, the VCS amplitude of electroproduction of real photons competes with a large and 
often dominant Bethe-Heitler (BH) amplitude, in which the real photon is emitted by leptons in the scattering 
process. In the leading order in electromagnetic interactions, single-spin asymmetries (SSA) in electroproduction of 
real photons are caused only by the VCS amplitude, including its interference with BH process. BH amplitude alone 
does not lead to SSA, unless higher-order electromagnetic corrections are included. Here we present a calculation of 
SSA coming from such corrections. 

The fact that QED loop radiative corrections can induce the beam SSA in 7 — p and — p collisions with production 
of a e + e _ -pair or radiation of a photon is known for a long time. A spin-momentum correlation in differential cross 
sections for these processes were first considered by Olsen and Maximon in 1959 0, and somewhat later in other 
publications [|| . It was shown that a large azimuthal asymmetry originates from interference of the first and second 
(two-photon exchange) Born approximation. These studies addressed primarily the case of low energies, such that 
lo/M << 1, where u> is the energy of the incoming polarized particle and M is the target mass. 

More recently it was pointed out, for example, in Ref.Q that SSA can be induced also by a pure loop correction to 
the lepton part of the interaction. For the related processes of radiative Moller scattering and electron-positron pair 
production on an electron target, the calculation of Ref.Q yields substantial SSA, reaching tens of per cent in selected 
kinematics. Asymmetries of this magnitude could present a significant systematic correction to VCS studies, and this 
fact partly motivated our present calculation spanning a broad kinematic range. In Ref.|8| the corresponding effect 
is included in the full radiative correction to beam SSA due to interference between BH and (the absorptive part of) 
nucleon VCS amplitudes. In Ref. 8] the full correction to the beam SSA was computed numerically, but up to now 
there is no published analytic expressions for the beam SSA caused by the loop correction to the leptonic tensor, that 
is simple enough to include in Monte Carlo generators for analysis of VCS experiments. The present paper aims at 
closing this gap. In addition, using an unpolarized leptonic Compton tensor originally derived in Ref.[9j, we obtain 
new results for SSA in the BH process with a polarized proton target arising from QED loop corrections. 

The QED effects considered in this paper are referred to as model- independent, since they do not require additional 
knowledge of the nucleon structure. They represent systematic corrections to SSA measurements such as |l0|, [ill, [l2|, 
0,0] in the interference region between BH and VCS amplitudes of electroproduction of real photons. 

II. GENERAL FORMALISM AND BEAM SINGLE-SPIN ASYMMETRY 



The contribution to the beam SSA for the BH process, 

e-(fci) + P(px) -» e-(fe) + 7 (fc) + P{P2), 



(1) 



2 




FIG. 1: One-loop corrected diagrams that induce single-spin asymmetries in the BH process (1). Only 
radiation from the outgoing electron contributes. Vertical solid lines indicate the unitarity cuts. 



induced by one-loop corrections to the leptonic part of interaction in the case of longitudinally polarized incoming 
electrons can be written in terms of contraction of the leptonic and hadronic tensors 



4tt Ba V H, 



.4" 
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where the symbol 5Re denotes the real part and B^ v is an unpolarized leptonic tensor for large-angle photon emission 
in the process (1). We define it as 



= ^Tr(k 2 + m)0A M (fci + m)O v \ , 

k 2 + k + m ki — k + m 
0\n = 7a 7^ + 7^ 7 7a , s = 2kk 2 , t = -2kk x 
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where m is the electron mass. Note that both Mandelstam invariants (s and t) in experiments that study proton VCS 
[Tfl ITU IT2I flit IT^j are large, therefore in our calculations we can omit the lepton mass in the quantity B^ v and write 
the latter in the form 

(s + u) 2 + (t + u) 2 4A 2 - - 
= — g^A —(k 1 f J ,kiu + k2 fl k 2v ) ,u = -2k x k 2 , A = k 2 -ki+k = Pl - P2) A 2 = s + t + u, (4) 



where the tilde notation for 4-vectors denotes the gauge-invariant substitution, a p 
For the hadronic tensor, we use its Born expression 



H MV = -Tr{p 2 + M)r M (p 1 + M)(l - l5 S)T„ , T M = (F 1 + F 2 ) ltl - P±t±p>± F 2 , 

where F\ 2 = Fi 2 {A 2 ) are the Dirac and Pauli proton form factors, respectively, and S is a 4-vector of proton 
polarization. The expressions for spin-independent and spin-dependent parts of the hadronic tensor are 



TS 
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where (fivAp) = e M „ApAA and the sign convention for the Levi-Civita tensor is £0123 = 1- Only the spin-independent 
hadronic tensor contributes to the beam SSA, therefore contraction of the tensors reads 



[s + u) 2 + (u + tf 
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where V = 2k\p\ : X = 2k 2 p\. 
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Let us now consider the leptonic tensor. In general, its spin-dependent part for the case of the beam longitudinal 
polarization can be written as 



P = P(°) + — P(i) (8) 



where the Born value PjS) is pure imaginary and antisymmetric and therefore it does not contribute in the beam 
SSA, whereas the one- loop correction to it Pjh) , together with the imaginary antisymmetric part, contains also the 
the real and symmetric part. The latter is caused by interference between the Born-level BH amplitude and the 
diagrams in Fig^with an additional photon loop and an external photon coupling to the outgoing electron. These 
diagrams produce a nonzero imaginary (absorptive) part of the amplitude in the physical region of the process (1) 
due to two-particle electron-photon intermediate states as shown in Fig^by the unitarity cuts. 

Calculation of the one-loop diagrams of Fig contributing to the radiative leptonic tensor can be done by standard 
QED techniques. Details of the calculation can be found in Ref.j^] for unpolarized electrons and in Ref.^ij for 
longitudinally polarized electrons. In calculations of SSA, we are only concerned with the absorptive parts of the 
electron virtual Compton amplitude that enters the radiative leptonic tensor. It should be noted that calculation of 
SSA is infra-red safe, i.e., infra-red singularities explicitly cancel at the loop level for this observable. It is an important 
consistency check of the calculation, since ultra-soft photon radiation that normally cancels such singularities has zero 
absorptive part and therefore cannot assist with infra-red divergence cancellation in the considered observable. 

The contribution to P^]) we are interested in can be derived by standard loop computation techniques of QED. 
Neglecting terms that are explicitly anti-symmetric in indices fx <-> is, we arrive at the following expression (c.f. 
Eq.(36) of Ref.[3), 

P$ = i(k 1 k 2 Av)[B 1 k lll + B 2 k 2 ^\ - i{k x k 2 An)[Blkx v + B* 2 k 2u ] , (9) 
where the quantity B\ has the form 
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where, except for the 4-vector A, we used the same notation as in ^(| (see also 0), 

(kik 2 Aa) = e a p\akiak 2 /3A\ , a = s + t, b = s + u, c = t + u, (11) 



- 2 A 2 , ^ t 



G = L Au (L A + L U - 2L t ) - y - 2Li 2 (l - ^-) + 2Li 2 (l - , 
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L xy = L x — L y , L A = In — — , L u = In — ^ , L t = In — % , L s = In — — 

m z m z m z m z 



The quantities B 2 and G can be derived from B\ and G by substitution 

B 2 = -B 1 (s^t), G = G(s++t). 

Imaginary parts of B\ and B 2 which induce the real symmetric part of P^J can arise from the terms containing 
L s and G, and the terms contributing to the imaginary part of G are L s and Li 2 (l — s/A 2 ) because of the condition 
1 — s/A 2 > 1. From the form of proparators in the Feynman diagrams of Fig^it follows that to obtain the imaginary 
part of B\ and B 2 , one has to add a small negative imaginary part to the electron mass. It leads to 

l-s/A 2 

s , s s f dx In x u + 1 
L s = \n— -nr, Li 2 {l-—)=- / \n——L sA , (12) 

o 

and it means that 

cj TO G = 27 rln^, 
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where the symbol Sm stands for the imaginary part. By combining the previous results, we arrive at 



2tt - ^ „ 2tt 
st st 



3mBi = Bx , SmB 2 = -B- 



2A 2 (A 2 -t) r3t + 2u 2, , t s , 

Bi = t - tww-t ln{1+ ^> (13) 

B 2 = 2 *.[tZ*i-2{l-»)Ml + ±)] 
t L t + u V t' K u n 

Note that the terms containing anomalous poles at t — ► are cancelled explicitly in Eqs. Moreover, both Bi 

and B 2 are proportional to t if t — > 0. Contraction of the tensors in the numerator of Eq.Q can be written in the 
following form that appears as a factor in the spin-dependent numerator for beam SSA: 

P«iV = ^^f^ (Ft - [(2V — s + A 2 )_Bi + [2X - s - u)B 2 ] . (14) 

An important observation is that the expressions l|13|l contain an overall factor of A 2 , that is a squared 4- momentum 
transferred to the proton target. Since experiments on DVCS select small values of A 2 , it directly implies additional 
suppression of beam SSA in this region. 

Equations ©,0, 113I14|1 determine the beam SSA for the BH process Q. Note that the expressions for SSA do 
not contain large logarithms involving the lepton mass. 

Let us now introduce kinematic invariants used to define experimentally measured asymmetries in VCS. The 
following set of kinematic variables |l5j is accepted for experimental analyses: Azimuthal angle $ between planes 
(q 1; ki) and (qi, P2) in laboratory system (with axes OZ along direction — qi = ka — ki), and invariant variables 

-q\ 2gigi 2 

y = — — , qi = fci - k 2 , q x = u. (15) 



2pm 1 " V 

The advantage of these variables is that one may use both the laboratory system and c.m.s. of subprocess 7(91) + 
P(pi) — > r y(k) + P(p 2 ) to investigate the azimuth correlation. The reason is that c.m.s. can be reached from laboratory 
system by a boost along the direction of qi, and such transformation does not change the azimuthal angle $. Therefore, 
we can use the simple expressions for the particle energies and their angles in c.m.s. 

pi = (Si, 0,0, |pi|), p 2 = {E 2 , |p 2 | sin6>Arcos$, |p 2 | sin0jvsin$, |p 2 | cos6»jv), (16) 
fci = £i(l,sin6'i,O,cos0i), k 2 = £ 2 (1, sin# 2 , 0, cos0 2 ), 
k = (k ,-p 2 ), qt = (ex - s 2 , -pi) 

in terms of variables l|15|) to form all the invariants that enter into beam and target asymmetries. They read 

V(l-xy) V(l-y + xy) V(y-xy) 

£ 1 = — n rrrr > £ 2 = TTIttt > ^0 = 



2y/sy 2y/ty 2-s/ty 

V(y + 2r) V(y-xy + 2r) 

fi l = — n rm — > -^2 



2y/£y 2V£V 

y(l-2x + xy + 2xr) y(l - 2x - y + xy - 2xt) 

cos6»i = — , cos6» 2 = — , (17) 

\{l-xy) \(l-y + xy) 

. a zVv&y . a 2\/v&y 

Smt>x = rry r , Sm^ 



A(l - xy) ' \{l-y + xy) 



y(y-xy + 2xt) - 2jp . 2£yjp+ - p)(p - P-) 

cos On — 77 r , sint/jy — 



X(y-xy) ' ' X(y-xy) 

where we used the following brief notation 

M 2 

A 2 = y 2 + Axyr, r]=l-y- xyr, £ = y - xy + t, t = — , (18) 
- p ± = ^-[(l-x)(y±\) + 2xr}. 



V ' r ^ 2C 
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Here we introduced a dimensionless variable p with the quantities p± having a meaning of the minimum (/0_) and 
maximum (p+) value of p at fixed x, y and V. By using relations l|17|) it is straightforward to obtain 

Vv 

s = — f-[2A'cos$ + xy(l + 2xt) + ^(1 - y - 2x + xy - 2xt)] , (19) 
Vv 

t = - — %-[2K cos$ + xy(l - y- 2xt) + p(l - 2x + xy + 2xt)] . 

/ 1 j a \ V 2 Ky . 2 ^a;(p + - p)(p - p_) 

(fcifc 2 Api) = — — sm$ , X = . 

2A y 

Our estimation of the beam SSA is demonstrated in Fig. 2 for the conditions of Jefferson Lab and HERMES 
experiments on beam SSA in VCS at different electron beam energies: Eb — 4.25 GeV [ljj, 5.75 GeV 0] an d 27.5 
GeV [ll| and fixed values of —qf = 1.25, 1.08 and 2.6 GeV 2 , respectively. It can be seen from Fig. 2 that the asymmetry 
is rather small, not exceeding 0.1 per cent in the kinematics of the considered experiments, even for a rather broad 
range of variable p. The calculated effect is smaller for the kinematics of planned SSA measurements in DVCS by 
COMPASS collaboration at CERN 0], since it covers smaller values of A 2 . The effect for kinematics an d 
are similar in magnitude. 




FIG. 2: Dependence of the beam SSA (in per cent) on azimuth angle (in degrees) for three different experi- 
mental conditions (left panel): CLAS1 10] corresponds to x = 0.19, y = 0.825, p = 0.024, for CLAS2 
x = 0.18, y = 0.5, p = 0.0185, and for HERMES 4l] x = 0.11 y = 0.458, p = 0.005. Curves for different 
values of p on the right panel correspond to CLAS2 conditions |12|l. 



III. TARGET SINGLE-SPIN ASYMMETRIES 



Let us now consider the target SSA in the process (JTJ caused by one-loop corrections to the unpolarized part of 
the leptonic tensor. In contrast with the beam single-spin correlation, where the effect is due to the real symmetrical 
part of the spin-dependent leptonic tensor, this time it is related with the imaginary antisymmetric part of the 
spin-independent leptonic tensor. In this case 

where the effect is induced by the imaginary and antisymmetric part of the spin-independent tensor K^J that can 
be computed directly (c.f. 0). Keeping only antisymmetric terms in p «-> v, the result can be written as 

K $ = ^ 5m [( T i 2 - T 2i)(K~ k 2» ~ hvhn)] , (21) 
where, neglecting terms proportional to to 2 , we have 



T 2 



A 2 , . A 2 , a9 ~ A ,/itA 2 2u-.s + t v 2us 

-5-c(it - s)G + ^-(uA 2 - st)G - 2 A 2 ( + —L Au ) + 8u + U - s + 

s z t z st a c 
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FIG. 3: The target SSA as a function of azimuthal angle For left panel, the polarization 4-vectors are 
defined by Eqs. 126H and the right panel corresponds to Eas. 1281 . 



4(w 2 - cs)(A 2 L Au - a) A 2 (2c + t)(st - uA 2 )L As A 2 c(2u- s)L At 



c 2 t 



bs 



(22) 



and T-ii is derived from T12 by substitution T21 — Ti2(t <-> s). 

After extracting absorptive parts from L s and G we obtain the following expression for antisymmetric imaginary 
part of the leptonic tensor l|21l) : 



st 



2st 



1. u + t 1' 



/1 u + t t\ 
4u -In 

\t u cJ 



(23) 



We note an overall factor of A 2 from Eq. l|23|l that leads to additional suppression of target SSA in DVCS kinematics. 
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Contraction of the antisymmetric tensors in expression (|20|l for target SSA reads 

aM{Fi + F 2 )TG S , 

where the quantity G s depends on the target-proton polarization 4-vector S and has the form 



H s K d) 



(24) 



G, = 



Ft 



AM 2 ' 



F 2 )(k 1 k 2 AS) 



2M 2 



{p 2 S)(k\k 2 Api 



(25) 



It follows from Eq. 1251) that in general the one-loop correction to the leptonic tensor in radiative process (1) 
generates the target SSA due to all three possible orientations of the target polarization. Here we consider two possible 
conventions for defining directions of target polarization. First, one can consider the case when the longitudinal 
polarization (L) in laboratory system is along the electron beam direction of ki , transverse polarization (T) lies in the 
plane (ki,k2) and the normal (N) one is along the normal vector (ki x k 2 ). The corresponding polarization 4-vectors 
can be expressed via the 4-momenta |17j 



S 



'Vt 



qN 



2(pk 1 k 2 p 1 ) 



S 



k-2^ - (1 - y - 2xyr)k lfi - xyp lfl 
y / Vxyrj 
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Then we have 



"si 
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V 3 xy 
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- xyFi + — (xy(l - y - 2xyr) +ty(l + 2xt) + p) 
( F i - £ F 2) Ml -V)- P(l - xy) - (2 - y)t] 



4F 2 (fc 1 fc 2 A Pl 
V 4 xyr 



(27) 



with the same proton form factors F\ y2 as in Eqs.(5,6). 

In another convention, one may choose directions to define the polarization 3-vector in the lab system. If the 
longitudinal direction is taken along p 2 , the transverse one in the plane (ki,p 2 ) and the normal is along 3-vector 
p 2 x ki, then the corresponding polarization 4-vector can be written as [lg 



oL 



-2tA m - 



^Vt P {At + p) 



Si - -HO^i , C = P {1 - xy + t) - T(xy-tf 



p(Ar + + (p - 2r(t - xy)A t _ t - p(2 - xy + i)p Xpi 

VvCpWTp) 
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In this case 



Lto„ 




Vt(At + p) 
P 

V((4t + p) 



(kik 2 A Pl )(Fi + F 2 ) 
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(29) 



The target SSA in the BH process (1) is shown in Fig. 3. Our calculations indicate that beam and target SSA 
generated by loop corrections to leptonic part of the interaction in the BH process are small and for the considered 
experimental conditions they do not exceed 0.1 per cent. The reason is that in addition to being multiplied by the 
fine structure constant a = y^, they contain additional suppression for the relevant values of kinematic invariants. 



IV. SUMMARY 



In conclusion, let us discuss the role of other radiative corrections in the BH process coming from real-photon 
radiation. In VCS experiments, the kinematic cuts are imposed in such a way that the phase space of the (undetected) 



8 



additional photon is restricted to its relatively small values. In this case, the main contribution to the radiative 
correction comes from spin-independent soft photon emission that does not affect polarization observables, but does 
change unpolarized cross sections by as much as about 20 per cent (see, e.g., @ for VCS case and Ref.hsLhgl for 
elastic electron-proton scattering). Therefore in such a soft-photon-cmission regime, the loop correction considered 
here is the only model-independent radiative correction to SSA. 

Thus, we demonstrated that systematic corrections to beam and target SSA arising from the higher-order QED 
effects in the BH process are negligible compared to the relatively large (tens of per cent) experimentally observed 
asymmetries 1 1 1| due to interference between the BH and VCS amplitudes in electroproduction of real photons. 
We confirm that the basic assumption of negligible SSA from the BH process alone holds to better than 0.1 per cent 
accuracy, thereby justifying present interpretation of SSA as arising mainly from the BH-VCS interference and VCS 
mechanisms. 
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